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We present results for the mass spectrum of charm and bottom heavy baryons, using MILC coarse 
lattice configurations with 2+1 flavors. Clover heavy quark propagators with the Fermilab inter- 
pretation and improved staggered light quark propagators are used to construct two point functions 
with local operators for different flavor and spin states. 
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1. Introduction 

Heavy baryons have been extensively investigated by experimental and theoretical approaches. 
From experiment, the singly charmed heavy baryon mass spectrum is well known; however, the 
other heavy baryon masses are only crudely known. Recently, from the DO [1] and CDF [2] collab- 
orations, there are measurements of the mass of the singly bottom baryon S^. From lattice QCD, 
there are several quenched calculations for the heavy baryon mass spectrum [3-8] , and most results 
are in fair agreement with observed values. In this project, we apply lattice QCD with dynamical 
sea quarks for the same objective. 

Previously, we have presented the preliminary results for the mass spectrum of singly charmed 
heavy baryons [9]. Here, we present the singly charmed and bottom heavy baryon as well as the 
doubly charmed and bottom heavy baryon mass spectrum. In this work, we use two different 
interpolating operators studied in Ref. [3], and construct two-point functions using the method of 
Wingate et al. [10] to combine staggered propagators for the light valence quarks and a Wilson 
type (clover) propagator for the heavy valence quark. In addition, we show that there is no taste 
mixing problem by studying the explicit spin and taste structure of operators. 

2. Construction of two-point functions 

The interpolating operators to describe singly heavy baryons are [3] 

^5 = Sabci^CYs^^H, = EabcWfCyMn, (2.1) 

where e a b c is the Levi-Civita tensor, yi and y/2 are light valence quark fields for up, down, or 
strange quarks, is the heavy valence quark field for the charm or the bottom quark, C is the 
charge conjugation matrix, and a, b, and c are color indices. Basically, G5 is the operator for 
s K = + , and is for s 71 = 1 + , where s K is the spin parity state for the light quark pair. As can 
be seen, the spinor indices of the light quark fields are contracted together, so the spinor index of 
the operator comes from Furthermore, each heavy baryon is obtained by choosing appropriate 
quark flavor combinations, and, for doubly heavy baryons, we can simply interchange the heavy 
quark field and the light quark fields in Eq. 2.1 [4]. 

Since we use staggered fermion propagators for the light quarks, and a Wilson type propagator 
for the heavy quark, we will use the method of Wingate et al. [10] to convert to a naive quark 
propagator starting from a staggered propagator. They show 

G Y (x,y) = a(x)tf(y)G x (x,y), (2.2) 

where 

n 

G v (x,y) is the naive quark propagator, and G x (x,y) is the staggered quark propagator. With this 
basic relationship, we can construct the heavy baryon two-point functions for s n = + using the 
operator ^5. 

C s (p,t) = Ee-^VsfrO^s&O)) (2-4) 
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= £ e-^e^e^MGf 7 (x, 0)Cy 5 Gf (x, 0)Cya]Gg (x, 0) (2.5) 
= £^ 1 > ? £ fl/x .£^ c 4G^(x,0)Gg(^,0)G| / c '(x,0). (2.6) 

X 

G\(x,0) and G2(x,0) are naive light quark propagators, Gjj(x,0) is the Wilson type heavy quark 
propagator, and G\ x (x,0) and G2%(x,0) wee staggered light quark propagators. Note that the trace 
in Eq. 2.5 is for spinor indices, not for color indices. Similarly, we can derive the two-point function 
for s n = 1 + with the operator 6^, which is 

C M v(/M) = £ e -^e a6c e flW 4(-l)^5^ v G^'(x,0)G^'(x,0)G^(x,0)^ (2.7) 

X 

We note that 8^ v in Eq. 2.7 is an artifact of using Eq. 2.2 to convert the light staggered propagtor 
to a naive propagator. With Wilson or clover light quarks, C^ v would have off diagonal elements. 
We shall now elucidate the origin of the 8 -function. 

3. Taste mixing 

We know that heavy-light mesons do not have any taste mixing problem [10]. Since we use 
two staggered quarks for singly heavy baryons, the taste structure of the operators might be more 
complicated than in the heavy-light meson case. Furthermore, presence of the 8 function in Eq. 2.7 
suggests we investigate more seriously. 

We have resolved this problem by considering the taste and spin structure of the operators 
explicitly. Let's look at the operator in Eq. 2.1. Since the heavy quark field carries no 
taste index, we can simplify the problem by considering the di-quark operator @„ rather than the 
original operator G^, i.e., 

% = wfC7n¥2- (3-D 
Next, we have a relation between the naive quark and the staggered quark fields 

¥ a '(x) = a a ' a (x) X a (x), (3.2) 

where \j/ a ' is the naive quark, % a are four copies of staggered quark, £2(x) is the matrix given in 
Eq. 2.3, a' is the spinor index of the naive quark, and a is the copy index. Also, we have 

X a {y + t;) = 2ti ia {t;)q ai - a {y) (3.3) 

and 

x = y + $, (3.4) 

where q m - a is the staggered quark in the taste basis, a is the staggered spin index, i is the taste 
index, x is the coordinate index for all lattice sites, y is the coordinate index for even sites only, and 
% is the vector of each corner of the unit hypercube. Using Eqs. 3.2 and 3.3, we can rewrite the 
naive quark y/ 01 (x) as a linear combination of the staggered quarks in the taste basis 

y a \x) = a a ' a ($)x a (y+Q = Cl a ' a 2a fia (£)q ai > a (y). (3.5) 
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Now, we can put Eq. 3.5 into the di-quark operator, and sum over the hypercube in order to get the 
continuum di-quark operator 3)™ nn {y), 

% onti (y) = KtfMcy^to) 

= £2n+ i «(^) g « i >(y)a w (§)(cr M )^'n^(^)2a^(^)^(^ (3.6) 

= £4Q t ' a (^)^( J )(-l)^(C7,) fl ^(^)/^(j), 

using the relation 

n rc '^)(Cy M r = (-l)^(Cy M y*. (3.7) 

Next we sum over % in Eq. 3.6, using 

£^' O! (^)(-l)^^(^)=4(C7 M ) a/3 0(7MC- 1 )' 7 . (3.8) 
5 

Finally, the continuum di-quark operator £F™ m '(;y) can be written 

%° nt \y) = l6q ai ^(y)(CY^ a P (Y.C-'yJq^iyXCY^. (3.9) 

From Ref. [11], we immediately notice that this continuum di-quark operator overlaps with 

a physical s 11 = 1 + spin parity state with a taste singlet. Thus, the original operator G^ overlaps 
with i and | total spin parity states of the heavy baryon without any taste mixing problem in 
the continuum limit. 

Similarly, we can easily derive the continuum di-quark operator for the operator G5 

gcon, i{y) = l6q au a{y){ c y5) ap {y 5 C- l )^q^- h {y){CY 5 ) a \ (3.10) 

and we can state that the operator G5 overlaps with physical | + total spin parity state of the singly 
heavy baryon with a taste singlet for the same reason. 

In addition, if we consider the two-point function C c ^ tl of the continuum di-quark operator 
@ c S n, \ then it appears 

C c ;f(y,0) = 16 2 TV[GiCy,0)(Cr M ) <8> (Cy M )+G 2 (j,0)(Cr v ) t <g> (Cy v )] 

X (CY^)ab{YvC^ l )b'a'8bb'Saa' (3.1 1) 

= 16 2 Tr[G 1 (y,0)(C7 / ,) (C7 M ) t G 2 (y,0)(C7 v ) t (C 7v )] 

xTr^CT,)^) 1 "], (3.12) 

where 

Tr[(C7 A( )(C7 v ) t ]=4<V- (3-13) 

As we can see from Eq. 3.12, C c ^""(y;0) consists of two traces. The first trace contains spinor and 
taste indices. The second trace is over copy indices. Now, we can address where the delta function 
in Eq. 2.7 came from: cancellations of the copy indices govern the delta function. Furthermore, we 
notice that even though the copy indices have a non-trivial structure, we can still use the operator 
G^, because the contribution of the copy indices is a constant, and independent of spin and taste. 
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4. Results 

We use three ensembles of 20 3 x 64 MILC coarse dynamical lattice gauge configurations with 
lattice spacing a fa 0.12 fm [12]. The details of the ensembles are shown in Table 1 . For each 
configuration, we require several propagators for the valence quarks. We compute nine different 
staggered light quark propagators with masses between 0.005 and 0.02, and three staggered strange 
quark propagators with masses 0.024, 0.03, and 0.0415, and each propagator is calculated with four 
different time sources. For the heavy quarks, we use hopping parameter values K = 0.122 for the 
charm quark, and K = 0.086 for the bottom quark based on tuning for our heavy-light meson decay 
constant calculation [13] [14]. 



J3 


ami 


am s 


# of conf . for charm 


# of conf. for bottom 


6.76 


0.007 


0.05 


545 


554 


6.76 


0.01 


0.05 


591 


590 


6.79 


0.02 


0.05 


459 


452 



Table 1: Parameters of the ensembles, and number of configurations used for each heavy quark. The light 
sea quark mass is m/ and the strange sea quark mass is m s . 

For the data analysis for the baryon propagators, we use the same fitting and extrapolation 
method as in our previous work [9]. 
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Figure 1: The mass splittings for singly charmed heavy baryons. 



We present mass splittings in MeV for singly charmed heavy baryons in Fig. 1 . Crosses corre- 
spond to the experimental data [15], and squares and diamonds with error bars indicate our results. 
The experimental errors are considerably smaller than our own, so we ignore them here, and the 
eiTors of our results are statistical errors. Since we do not have partially quenched staggered chiral 
perturbation theory for singly heavy baryons to fit both sea and valence quark mass dependences 
simultaneously, we use four slightly different extrapolation methods. The result from each method 
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is labeled as Full, Full2, Partial, or Partial2 in Fig. 1. The first method takes all partially quenched 
data, and linearly extrapolates for the valence quark masses and then the sea quark masses; this 
is the "Partial" result. Alternatively, we consider only data where light valence masses equal the 
corresponding sea quark mass and do a single extrapolation; this is the "Full" data. Further, we 
can obtain the mass splittings by a second method. Instead of applying the chiral exptrapolations 
to the masses themselves, we apply them to the splittings. Results from the second procedure are 
denoted with a "2". 
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Figure 2: The mass splittings for singly bottom heavy baryons. 



We also present mass splittings for singly bottom heavy baryons in Fig. 2. Crosses are the 
recent experimental data from the CDF [2], and octagons are the experiment data from the D0 [1]. 
We include error bars for the D0 results, since the errors of the D0 results are larger than the other 
experiments, about 1 8 MeV. 
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Figure 3: The mass spectrum for doubly charm (a), and bottom (b) heavy baryons. 
In addition, we present a preliminary mass spectrum in lattice units for doubly bottom and 
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charmed heavy baryons in Fig. 3. 
5. Future study 

We have presented the mass spectrum for charm and bottom heavy baryons based on dynam- 
ical lattice QCD, and a solution of the taste mixing problem of the operators. The study can be 
extended in several ways. In the near future, we would like to extend this work into 28 3 x 96 MILC 
fine lattice gauge configurations; so that we can determine the mass spectrum more accurately. In 
this work, we apply the simple linear extrapolations for sea quark masses and valence quark masses 
separately; as a result, our partially quenched data analysis gives larger errors than the errors from 
full QCD data analysis. In order to extrapolate for the both sea and valence quark masses simulta- 
neously and accurately, we would like to develop partially quenched staggered chiral perturbation 
theory for heavy baryons. Next, we must consider systematic errors, such as finite size effects and 
discretization errors. In addition, we are interested in the | excited state masses of singly heavy 
baryons, which cannot be studied with our current method. 
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